Measurement of the electronic compressibility of bilayer graphene 
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We present measurements of the electronic compressibility, K, of bilayer graphene in both zero and finite 
magnetic fields up to 14 T, and as a function of both the carrier density and electric field perpendicular to the 
graphene sheet. The low energy hyperbolic band structure of bilayer graphene is clearly revealed in the data, 
as well as a sizable asymmetry between the conduction and valence bands. A sharp increase in A" -1 near 
zero density is observed for increasing electric field strength, signaling the controlled opening of a gap between 
these bands. At high magnetic fields, broad Landau level (LL) oscillations are observed, directly revealing the 
doubled degeneracy of the lowest LL and allowing for a determination of the disorder broadening of the levels. 
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The hyperbolic band structure predicted for bilayer 
graphene stands apart from the quadratic dispersions com- 
mon to almost all semiconductor systems at low energies, and 
is equally distinguished from the linear dispersion of mono- 
layer graphene [1]. The curvature near the band edges rapidly 
evolves toward a near-linear form with increasing energy [2]. 
Accordingly, the single particle density of states (DOS) has 
an unusual form, rising linearly with increasing energy from a 
finite value at zero energy. Owing to the unique, tunable band 
gap of bilayer graphene, the DOS can be sharply modified re- 
flecting the changing band structure [3, 4]. Here we have used 
an electric field penetration technique to measure the com- 
pressibility, K, of the two-dimensional electronic system in 
bilayer graphene, and thereby observe and explore changes 
in the DOS as a function of applied electric and magnetic 
fields [5]. Thermodynamic quantities such as the compress- 
ibility are sensitive to both localized and extended electronic 
states, and are known to be an effective tool for investigat- 
ing many-body effects that arise in interacting systems and 
have recently been predicted to occur in bilayer graphene [5- 
8] . Our experimental results demonstrate that, contrary to the 
usual parabolic approximation made at low energies, the hy- 
perbolic nature of the band structure must be accounted for to 
the lowest carrier densities measured. Moreover, a clear asym- 
metry between the conduction and valence bands is revealed, 
implying a reduced DOS in the conduction band relative to the 
valence band. When a strong electric field is applied perpen- 
dicular to the graphene sheet, a rising peak at zero density is 
consistent with the opening of a band gap. At high magnetic 
fields, B, strong oscillations in the compressibility signal arise 
due to the formation of Landau levels (LLs), which are con- 
siderably broadened by disorder. 

Bilayer graphene is composed of two monolayer graphene 
sheets in AB (Bernal) stacking. To first order in the tight- 
binding (TB) approximation, the electronic dispersion of pris- 
tine, ungated bilayers consists of two nested hyperbolae in 
each of the conduction and valence bands, occurring at both 
the K and K' points of the Brillouin zone [1, 2]. The two 
lowest energy bands meet at k = 0, while the two higher 
bands are split off by a large interlayer hopping energy of 
about 7! = 0.4 eV. This simple yet highly unusual disper- 



sion can be radically altered by the application of an electric 
field normal to the plane, which breaks inversion symmetry 
and results in the opening of a tunable band gap between the 
lowest bands [3, 4, 9-12]. In a dual-gated field-effect transis- 
tor device, where a bilayer is sandwiched between two electri- 
cal gates, the Fermi level (and hence carrier density) and the 
size of the gap can be independently manipulated [9, 10]. Be- 
sides holding promise for future applications, this allows for 
investigation of physics depending on only one or the other 
quantity. With the bilayer embedded in a dielectric medium, 
the density and gap are controlled by the back and top gate bi- 



(a) p(kQ/D) 

20 40 60 



(b) 




-2 

-40-20 20 40 60 

V b (V) 




(c) 



top gate 

dt 

graphene : - 
back gate ^ - 



30 \irc\ 



-II— 

It oc Cm 
'V, 



v h 4 



FIG. 1: (Color online) (a) Sheet resistivity of the bilayer graphene 
device vs. top and back gate voltages Vt and Vb. The arrows indicate 
the orthogonal axes of the carrier density, n, and the displacement 
electric field perpendicular to the sheet, D av . (b) Optical microscope 
image of the device, (c) Schematic of the electric field penetration 
measurement. 
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ases, Vt and V b - Electrical gating of the graphene results in a 
carrier density n = a(AVb — f3AV t ), while the average value 
of the electric displacement field perpendicular to the sheet, 
D av = (D b + D t )/2 = (e b AV b /d b -e t AV t /d t )/2, generates 
the band gap. Here AVi = Vj — Vm, where Vm is an offset 
voltage from charge neutrality (CN) due to extrinsic dopants, 
while Ci and di are the dielectric constant and thickness of 
each gate insulation. 

Our dual-gated devices begin with Scotch tape exfoliation 
of graphene bilayers onto Si/Si02 wafers having a 310 nm 
thick oxide. Standard electron beam lithography and thin film 
deposition techniques are used to fabricate several Cr/Au con- 
tacts to the bilayer. A layer of dilute 950K PMMA is then 
deposited and overexposed with an electron beam dosage of 
24,000 ^C/cm 2 to create a thin (d t = 25 - 30 nm), hard- 
ened dielectric layer of amorphous carbon over the graphene 
[13, 14]. In the final fabrication step, a Cr/Au top gate is de- 
posited, closely aligned to the contours of the graphene. An 
optical microscope image of a typical finished device is shown 
in Fig. 1 (b). The mobility of the sample used in this work is 
2,900 cm 2 /Vs, determined from the linear density dependence 
of the conductivity along the line D av = [15]. All measure- 
ments were carried out at T=1.5 K. 

Electronic transport measurements were made using stan- 
dard techniques and serve to identify the bilayer nature of the 
sample, as well as to calibrate a and f3. In Fig. 1 (a) we 
show a map of the sheet resistivity at zero magnetic field as 
a function of V b and V t . In this plot, the charge carriers are 
electrons to the right of the diagonal and holes to the left. The 
peak resistivity occurs along the diagonal corresponding to 
CN; note the offset such that for V b = Vt = 0, the bilayer is 
hole-doped. Along this CN diagonal, the minimum resistivity 
occurs at the saddle point located at (Ho, Vto) = (—0.5, 1.9), 
which we identify as where both n and D av , and hence the 
induced band gap, are zero. Moving along the diagonal away 
from this point, the resistivity increases monotonically due to 
the opening of the band gap. Together with the anomalous 
quantum Hall effect peculiar to bilayer graphene and observed 
in this sample, these transport measurements demonstrate that 
our sample is a graphene bilayer [2, 16]. We establish the gate 
voltage to density calibration, a = 7.1 x 10 10 V _1 cm -2 , by 
counting oscillations in the compressibility data with chang- 
ing LL filling factor, v, at finite magnetic field, as in Fig. 3 
(a). Here v — nh/(geB), where g = 4 is the combined spin 
and valley degeneracies. The slope of the CN diagonal in Fig. 
1 (a) gives the ratio of the back to top gate capacitances per 
area, 1//3 = —d t e b /(d b e t ) = —0.11. Thus for SiC>2 with 
e b = 3.9eo we find e t = 3.3eo for the dielectric constant of 
the overexposed PMMA. 

The compressibility is measured using an electric field pen- 
etration technique schematically depicted in Fig. 1 (c) [5]. 
In practice, we measure the capacitance of the back gate to 
the top gate, C bt , where the electric field lines between the 
gates penetrate the grounded graphene sheet. The finite DOS 
limits the ability of the graphene to screen the electric field 
lines, establishing a direct relationship between the measured 



quantity, C b t, and the inverse DOS, d/i/dn, which is related 
to the electronic compressibility via K~ x = n 2 dfi/dn where 
(i is the chemical potential. For the device configuration used 
here, in the limit of an infinitely thin graphene sheet the mea- 
sured capacitance is given by 

r A t d g e b e t 

a g {etd b + e b d t ) + e d b d t 

Here A t is the area of the top gate and d g = (e /e 2 )<9/i/<9n, 
with — e the electron charge, is a convenient parameterization 
of dfi I dn in units of length. For typical experimental parame- 
ters, d g s=s 2 A. Referring again to Fig. 1 (c), with the graphene 
held at ground we apply a small AC signal (typically 0.5 V 
at 100 Hz) to the back gate and record the resulting current, 
I t oc Cbt, flowing from the top gate. Thus the quantity we 
are sensitive to is dfi/dn. The parallel stray capacitance, C s , 
arises from fringing fields both at the edge of the graphene 
flake, and between the top gate wiring and the entire Si sub- 
strate. Shielding the top gate wiring with an on-chip ground 
plane, visible in Fig. 1 (b), and soldering the top gate wire 
to a small diameter coax cable held less than 1 mm above the 
device help to minimize the impact of C s . The quadrature sig- 
nal in all measurements is small and constant across a range 
of frequencies / = 30 — 1000 Hz, confirming the capacitive 
nature of the measurement. 

Figure 2 (a) displays a map of the measured C bt at zero 
magnetic field, shown as a function of both Vt and V b . Mim- 
icking the sheet resistivity, the signal is strongest along the 
CN diagonal and monotonically declines with increasing hole 
or electron density. For clarity, Fig. 2 (b) and (c) show two 
profile cuts through the data, taken along the lines of constant 
D av marked in Fig. 2 (a) and plotted as a function of the 
carrier density. These mountain-like profiles evince a prime 
result of this work: since Cu oc dfi/dn, if any portion of 
the band structure were parabolic, a constant C bt vs. n would 
result. Such a parabolic approximation is often assumed for 
carrier densities <4x 10 12 cm~ 2 . Yet in each trace, the data 
monotonically and sharply decrease from a peak value near 
zero density. 

Insight into the decrease of dfi/dn with \n\ can be gained 
from a minimal tight-binding model of bilayer graphene. In 
the most basic TB model, the low energy bands of ungapped 
graphene give rise to a d/i/dn that peaks at n = and falls 
off roughly as 1/n [1, 2]. Utilizing this model with the Fermi 
velocity set to vp = 10 6 m/s, along with Eq. 1 where C s is 
the only fitting parameter, we fit the D av = trace with the 
dotted line in Fig. 2 (b). We take C s = 54.8 fF, such that the 
remaining signal varies by roughly a factor of two over the 
range of densities measured. Remarkably, this simple model 
reflects the general trend of the data, strongly suggesting the 
low energy band structure is in fact hyperbolic. While an- 
ticipated theoretically, this low energy regime has not been 
clearly observed previously. 

A clear electron-hole asymmetry with higher Cu values for 
electron doping is visible in both Fig. 2 (b) and (c), which is 
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FIG. 2: (Color online) (a) Measured back gate to top gate capacitance Cm, proportional to d/i/dn, plotted vs. Vt and 14 at zero magnetic 
field. The solid and dashed lines are profile cuts through the data while the insets schematically depict the band structure, for (b) D av — 
(zero gap) and (c) D av = 0.4 V/nm (gapped), respectively. The dotted, dashed, and solid lines in (b) and (c) are fits to the data described in 
the text. 



absent from the transport data. Both theory [17-19] and ex- 
periment [11, 12, 20, 21] expect that particle-hole symmetry 
is broken in bilayers leading to an asymmetry in the DOS. In- 
deed, in the full TB model several terms can introduce such 
an asymmetry [22]; the most straightforward method is to 
include a hopping velocity, V4, representing same-sublattice 
hopping between the two graphene sheets. Inclusion of this 
term increases (decreases) the curvature of the conduction (va- 
lence) band, thereby reducing (increasing) the DOS so that 
in the absence of disorder, du/dn undergoes a discontinuous 
jump between the two bands. Repeating the fit with a value 
of V4 = (0.6 ± 0.2) x 10 5 m/s yields the best result shown as 
the dashed curve in Fig. 2 (b). This V4 is close to the range 
0.4 — 0.5 x 10 5 m/s determined from spectroscopic experi- 
ments [12, 18, 20]. 

Generic to both mono- and bilayer graphene samples on a 
Si02 substrate are the so-called electron-hole "puddles" in- 
duced by charged impurities, which lead to inhomogeneous 
variations in the value of n across the sample over a typical 
length scale of tens of nm [23, 24]. Such fluctuations are av- 
eraged over in our measurement due to the large area of the 
device, so we include this effect in our model by convolving 
the fits with a Gaussian and adjusting the variance, Sn, until 
a best fit is found. The solid line in Fig. 2 (b) is found for 
Sn — 3 x 10 11 cm -2 . This Sn is consistent with other values 
in the literature for graphene-on-Si02 samples [23-26]. 

Overall, the fit is best at low densities, but becomes in- 
creasingly poor for higher hole densities suggesting our sim- 
ple model is incomplete. By including additional TB terms an 
improved fit may be achieved, but only at the cost of poorly 
constraining the fitting parameters. While the relatively low 
mobility of our device implies that effects beyond the single- 
particle band structure are secondary, it is likely that du/dn is 
impacted by screening and other many-body interactions not 
included in our model at this time. In fact, recent calcula- 
tions of the bilayer compressibility suggest interaction effects 
are not negligible [6-8]. In Ref. [6], a negative divergence 



of du/dn akin to that in parabolic 2D systems is predicted 
at very low densities. Although we do not observe this phe- 
nomenon, the calculation does not account for the presence of 
disorder. Alternatively, a recent calculation of interactions in 
the random phase approximation finds that du/dn increases 
monotonically with decreasing density, and is always larger 
than the non-interacting value [8]. While neither disorder nor 
an asymmetry effect are accounted for, the predicted increase 
in d^i/dn may qualitatively account for the greater variation 
seen in the data compared to the TB model fits. 

The profile cut for D av = 0.4 V/nm shown in Fig. 2 (c) 
exhibits a sharp peak centered at zero density, which we inter- 
pret as the opening of a band gap with concomitant changes 
in the DOS. Extending the minimal model used above to in- 
clude a gap suffices to describe this data trace. We first solve 
the Hamiltonian for a gapped bilayer and find the associated 
d[i/dn, accounting for V4 along with a band gap of size A. 
Then Cm is found using Eq. 1, and employing the same Gaus- 
sian smoothing over Sn as above we vary A until a best fit is 
obtained. The dashed and solid lines in Fig. 2 (c) are the fits 
before and after smoothing. For D av — 0.4 V/nm, the best 
fit is found with A = 26 meV. This is considerably smaller 
than the 140 meV "bare" band gap associated with a field of 
this strength applied between the sheets of the bilayer. Charge 
flow between the two layers in response to the external field 
will screen the bare gap [3, 19, 27, 28]. We note a gap of size 
A = 26 meV agrees to within uncertainty with that found 
in similar dual-gated devices at an equivalent displacement 
field D av as determined by IR spectroscopy [10]. Besides the 
density variations Sn, the DOS will in general be broadened 
by disorder or trigonal warping which will diminish the band 
gap further and smooth over sharp features [18, 29, 30]. In 
the fitting so far we have ignored these additional sources of 
broadening, so that the value of Sn used to obtain good fits to 
the data should be taken as an upper limit. 

In strong magnetic fields the compressibility is greatly 
modified by Landau quantization. In Fig. 3 (a) we display 



4 




FIG. 3: Color online) (a) Back gate to top gate capacitance Cm as a 
function of the top and back gate voltages, for B = 14 T. (b) Profile 
cuts along the line D av = for B — 0, 10, and 14 T, offset for 
clarity. Fits to the 14 T trace are for: density fluctuations only, 8n = 
3 x 10 11 cm -2 (dashed line); and Sn plus additional broadening of 
the LLs of width 8E — 11 meV (solid line). 



a map of Cm for a field B = 14 T, where the quantiza- 
tion alternately manifests as valleys, indicating compressible 
regions where a Landau level is partially filled, and ridges, 
when the Fermi level lies between a pair of LLs and the sys- 
tem is incompressible. The central valley is twice as wide 
as the others due to the unusual bilayer LL structure which 
places two 4-fold degenerate LLs at zero energy [2]. At 14 
T, the gap opening has only a weak effect on the LL struc- 
ture, apparent as an increase in Cm to the lower right of Fig. 
3 (a), with an accompanying decrease in the peak heights of 
the incompressible ridges corresponding to LL filling factors 
v = ±4. The asymmetry observed at B = T persists 
to high fields as stronger oscillations of Cm for electron vs. 
hole densities. Qualitatively, the increased curvature of the 
conduction band results in greater LL spacings which, for a 
given LL broadening, should result in a smaller remnant DOS 
(and thus larger Cm) in the gaps between LLs as compared 
to the valence band. In Fig. 3 (b) cuts through the data for 
B = 0, 10, and 14 T are plotted. The evolution of the asym- 
metry and the increased range of Cm values with increasing 
B are clearly visible. We construct fits to the data at 14 T 
by first assuming the underlying LLs are ^-functions centered 
at energies calculated within the minimal tight-binding pic- 
ture plus a term accounting for the V4 asymmetry [31, 32]. 
From this model DOS, Cm is calculated for the same values 
of C s and Sn as above. The resulting fit, shown as the dashed 
line in Fig. 3 (b), greatly overestimates the amplitude of the 
Cm oscillations. To improve on this, we have either to as- 
sume that Sn is nearly twice as large as the zero field value, 
or that an additional source of broadening is present. While 
Sn will not vary greatly with B field, disorder-broadening 
of LLs is expected in 2D systems including graphene [33- 
35]. The solid line fit in Fig. 3 (b) represents a calculation 
beginning with Lorentzian-broadened LLs all having equal 



width, SE =11 meV. This corresponds to a scattering time 
r = 60 fs, in good agreement with the transport scattering 
time r M = m*/i/e = 58 fs, where m* — ji/2vp. After 
further averaging over the same Sn as at zero field, a better 
agreement with the data is found. Broadening with Gaussian- 
shaped LLs consistently produces a worse fit. The asymmetry 
at high fields is not as well described by the zero field value of 
V4. As with the zero field fitting in Fig. 2 (b), the fits at high 
field match the data best at low carrier densities, and exhibit a 
similar discrepancy at high hole densities where the data fall 
below the fit. 

The picture emerging from our data is clearly consistent 
with a hyperbolic band structure, across a wide density range 
relevant for many experiments on bilayer graphene. Addi- 
tionally, a sizeable conduction-valence band asymmetry is 
present. Beyond the smoothing due to averaging over density 
fluctuations, in high magnetic fields additional broadening of 
the LLs is required to model our data. 

Note added in proof. Recently, we became aware of similar 
measurements of the compressibility of bilayer graphene by 
Young et al. [36] 
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